A GENERAL APPROXIMATION OF QUANTUM GRAPH VERTEX 
COUPLINGS BY SCALED SCHRODINGER OPERATORS ON THIN 

BRANCHED MANIFOLDS 



PAVEL EXNER AND OLAF POST 



Abstract. We demonstrate that any self-adjoint coupling in a quantum graph vertex 
can be approximated by a family of magnetic Schrodinger operators on a tubular net- 
work built over the graph. If such a manifold has a boundary, Neumann conditions are 
imposed at it. The procedure involves a local change of graph topology in the vicin- 
ity of the vertex; the approximation scheme constructed on the graph is subsequently 
'lifted' to the manifold. For the corresponding operator a norm-resolvent convergence 
is proved, with the natural identification map, as the tube diameters tend to zero. 
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1. Introduction 



The concept of quantum graph [EKK+08] serves as a laboratory to study quantum 
dynamics in situations when the configuration space has a complicated topology. At the 
same time, it is a useful tool in modelling numerous physical phenomena. To employ its 
full power, one should be able to understand the meaning of parameters associated with 
vertex coupling in such models, because one can typically associate many self-adjoint 
Hamiltonians with the same graph. An old and natural idea was to select plausible 
ones with the help of "fat-graph" approximations; it was formulated for the first time 
by Ruedenberg and Scherr [RuS53j who proposed a heuristic Green-formula argument 
to demonstrate that such a shrinking limit would yield the simplest coupling conditions 
conventionally labelled as Kirchhoff. After the interest to the problem had been renewed 
about twenty years ago a lot of effort was made to establish this limit rigorously, both 
for one-body Schrddinger e quation [FW931 ISaOOl iRgOTl IKuZOlL IEP051 IP06] and the 
(y-^ ' Ginzburg-Landau dynamics |RS01] . The idea of fat graph approximations has also been 
used in spectral geometry by |CdV86j in order to show that the first non-vanishing 
eigenvalue of a compact manifold of dimension three or higher can have arbitrarily high 
multiplicity. Moreover, (rescaled) fat graphs and their limits can be used in calculating 
spectral invariants — see [Gri08b| for a survey and [MM06J for an example of a graph 
with one edge and two vertices. We refer to the monograph [PI 2] for a detailed discussion 
of these problems and an extensive bibliography. 

However, already in the seventies investigations of branched electromagnetic wave- 
guides [Me78] indicated that the low frequency behavior, closely linked to the shrinking 
limit, can be different for different geometries, and after geometrically induced states 
bound states in Dirichlet tubes were discovered it became clear that the answer depends 
substantially on what boundary conditions one chooses for operators on the tube-like 
manifolds constructed ove r the g raph "skeleton" , and that the limiting coupling may not 
be of the Kirchhoff type ES89J. More recent investigations investigations have shown 
that this is typically the case for Dirichlet tube networks, cf. [P05l IMV071 Cr08l IACF071 
CE071 IDClOj . Here an energy renormalisation is needed and when one chooses the 
natural one which consists of subtracting the lowest transverse eigenvalue which blows 
up when the tube diameter e tends to zero, a nontrivial limit is achieved provided the 
fat graph from which one starts has a threshold resonance. 
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It is natural to ask whether and how can a shrinking limit produce all the admissible 
vertex couplings. Far from being a mere mathematical conundrum the problem is of 
practical importance. When building network-type objects one is primarily interested in 
control of the transport; the possibility to manipulate the junction dynamics is one of the 
most direct ways to achieve this goal. Methods to alternate transport properties of beam 
splitters by changing the junction geometry |CCL06[ IZLWllj or material characteristics 
|TY03j have been devised for photonic crystals, and it is only a matter of time when the 
fabrication technique progress will allow to address analogous technological challenges 
for networks of metallic nanowires or carbon nanotubes, modeled by by manifolds with 
Neumann boundary conditions or without a boundary, respectively, or semiconductor 
ones for which Dirichlet conditions are used. 

Furthermore, some applications of junction control can be made more specific. For 
instance, Cheon et al. |CTF04] proposed the generalized point interaction, i.e. the sim- 
plest nontrivial graph with two edges and a general vertex coupling, as a model gate for 
quantum computing; in a similar way other star graphs can model an arbitrary audit. 
The point is that the geometry of the eigenvalue manifolds of the corresponding opera- 
tors is described by the group U(n) where n is the number of the connected edges; for 
a concrete way how the corresponding eigenvalue anholonomy can be used in Grover 
search algorithm see, e.g. |TM07t ITN10] . To implement such proposals with real-world 
objects, a proper understanding of the junction dynamics is again essential. 

The aim of the present paper is to provide a complete solution to the problem for 
tubular network manifolds the boundary of which is either Neumann or absent. The 
approximation we are going to construct has several ingredients. The first is the use of 
scaled potentials. If the network dynamics is described by the Laplace-Beltrami operator 
the limit leads to the Kirchhoff coupling, hence one has to replace it by a suitable family 
of Schrodinger operators. One proceeds at that in two steps, first an approximation is 
constructed on the graph itself and subsequently it is "lifted" to the tubular manifold. 
In this way we have been able in [EP09j to approximate two important coupling types 
usually referred to as 5 and S' s . Referring to the graph approximation result obtained in 
[ET07] we conjectured existence of such approximation to any vertex coupling with real 
coefficients which covers all the couplings invariant with respect to the time reversal. 
We are going to show here that one is not only able to prove the said conjecture but in 
fact can do better: following the "algebraic" work done in |CET10j we demonstrate here 
existence of a "fat-graph" approximation for all self-adjoint vertex couplings. The new 
idea here is to lift the indicated, quite subtle approximation from the metric graph level 
to fat graphs keeping a precise control of the estimates. 

Let us recall briefly how the approximation constructed in [CETlOj works, a detailed 
description will be given in Section [2] below. It has several steps: 

(i) we change locally the graph topology disconnecting the edges and connecting 
the loose ends by addition finite edges the length of which tends to zero. Some 
of them may be missing, depending on the coupling we want to approximate 

(ii) the additional edges will be coupled to the original ones by 5 conditions of the 
strength dependent on the approximation parameter. We also add a parameter- 
dependent 5 interaction to the centre of these finite edges 

(iii) in order to accommodate the couplings with non-real coefficients we add mag- 
netic fields described by (the tangent components of) appropriate vector poten- 
tials, also dependent on the approximation parameter 

The main result of this paper consists of "lifting" this approximation to tubular net- 
works and demonstrating that one can approximate in this way any self-adjoint vertex 
coupling. Since the approximation bears a local character we concentrate our attention 
on star graphs having a single vertex; an extension to general graphs satisfying suitable 
uniformity conditions can be performed in the same way as in [P12] . 
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Let us remark in addition that in contrast to the approximation of Kirchhoff coupling 
by Neumann Laplacians on a tubular network the limit constructed in this paper is non- 
generic, cf. [ENZ01J, and at the same time non-unique. An example of different limits 
for the same coupling will be mentioned in Section I5.2[ and one can conjecture also ex- 
istence of significantly different approximation schemes, in particular, purely geometric 
ones |KP12j . The result of the present paper thus allows us to make several conclusions. 
From the physics point of view it answers affirmatively the question whether one can 
approximate all the vertex couplings allowed by the sole requirement of probability cur- 
rent conservation, and at the same time, it suggests one possible construction to achieve 
this goal technically. On the mathematics side, approximations of different vertex cou- 
plings open interesting possibilities in connection with the mentioned use of fat graphs 
in calculating spectral invariants. 

The paper is organised as follows: In the next section, we outline the approxima- 
tion procedure on the graph level. In Section [21 we construct the graph-like manifold 
model. Moreover, we introduce the quadratic forms corresponding to our operators on 
the graph and the manifolds and relate them with the "free" operators, i.e the corre- 
sponding Laplacians. In Section H] we briefly recall the convergence of operators and 
forms acting in different Hilbert spaces, apply the abstract conclusions to our situation 
here, and demonstrate our main result expressed in Theorem 14 .71 In Section [51 we present 
some examples, including the case of a metric graph embedded in M. u when the manifold 
model is an ^-neighbourhood of the graph. 

2. Approximation on the graph level 

As we have indicated in the introduction the approximation is constructed in two 
steps. First we solve the problem on the graph level, and the obtained approximation is 
then "lifted" to network-type manifolds. The first part of this programme was realised 
in |CET10j and we summarise here the results as a necessary preliminary. 

Any self-adjoint coupling in a vertex of degree n can be expressed through vertex con- 
ditions — one usually speaks about admissible conditions — which involve the boundary 
values /(0), /'(0) G C n . They are conventionally written in the form 

4f(0) + £/'(0) = 0, (2.1) 

where A, B are n x n matrices such that the n x In matrix (A\B) has maximum rank 
and AB* is Hermitian, cf. |KS99] . A pair (A, B) describing a given coupling is naturally 



not unique and there are various ways how to remove the non-uniqueness, see e.g. [Ha00[ 
IKu04] . The most suitable for our purpose is the one given by the following claim proved 
in |CET10j ; it is simple but it requires an appropriate graph edge numbering. 

Proposition 2.1. For a quantum graph vertex of degree n, the following is valid: 

(a) If S G C mxm with m < n is a Hermitian matrix andT G C mx ( n_m ) ; then the equation 

^ o) /,(0)= (-T* l(n-m))f(0) (2-2) 

expresses admissible vertex conditions which make the graph Laplacian a self-adjoint op- 
erator. 

(b) Conversely, for any self-adjoint vertex coupling there is a number m < n and a num- 
bering of edges such that the coupling is described by the conditions (I2.2p with uniquely 
given matrices T G C mx ( n - m ) and S = S* G C mxm . If the edge numbering is given one 
can bring the coupling into the form (12. 2 p by a permutation (1, . . . , n) (11(1), . . . , IT(n)) 
of the edge indices with the matrices S,T uniquely determined by the permutation II. 

Now we can describe the approximation of such a general vertex coupling. For sim- 
plicity we consider a star graph of n semi-infinite edges; in view of the proposition we 
may suppose that the wave functions are coupled according to (12 .2p renaming the edges 
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if necessary. The construction has two main ingredients. First of all, we have to change 
locally the graph topology, adding vertices to the graph as well as new edges which would 
shrink to zero in the limit. In this way one is able to get (12. ip with real matrices A, B; 
to overcome this restriction we need to introduce also local magnetic fields, i.e. to place 
suitable vector potentials at the added edges. 

The construction is sketched in Figure (TJ we disconnect the edges of the star graph 
and connect their loose endpoints by line segments supporting appropriate operators 
according to the following rules: 

k 




Figure 1. The approximation scheme for a vertex of degree n = 3 and 
n = 5. The inner edges are of length 2d, some may be missing depending 
on the choice of the matrices S and T. The arrows symbolise the vector 
potential. 



(i) As a convention, the rows of the matrix T are indexed from 1 to m, while the 
columns are indexed from m + 1 to n. For the sake of brevity, we use in this 
section the symbol n := {1, . . . , n}. 

(ii) The external semi-infinite edges of the approximating graph, each parametrised 
by s G M + are at their endpoints Vj connected to the inner edges by 5 coupling 
with the parameter Wj(d) for each j G n (see below). 

(iii) Certain pairs Vj, Vk of external edge endpoints will be connected by segments (or 
inner edges, labelled by {j, k}) of length 2d. This will be the case if one of the 
following conditions is satisfied, taking into account the convention (jl|): 

(a) j G rh, k > m+ 1, and Tjk ^ (or j > m + 1, k G rh, and Tjy ^ 0), 

(b) j,kem and (3/ > m + l)(T i; ^ A T kl ^ 0), 

(c) j, k G rh, Sjk 7^ 0, and the previous condition is not satisfied. 

(iv) We denote the centre of such a connecting segment by v {j^} and place there 
S interaction with a parameter w^j^jid). We adopt another convention: the 
connecting edges will be regarded as union of two line segments of the length d, 
with the variable running from zero at W{j t k} to d at Vj or Vk- 

(v) Finally, we put a vector potential on each connecting segment. What matters is 
its component tangential to the edge; we suppose it is constant along the edge 
and denote its value between the points and Vj as A^ k){d), and between the 
points V{j t k} and v k as A( k ,j)(d); recall that the two half-segments have opposite 
orientation, thus A^,j)(d) — ~^-(j.k){d) holds for any pair {j, k}. 

The choice of the dependence of Wj(d), W{j^}(d), and A^ k ){d) on the length parameter 
d is naturally crucial; we will specify it below. We denote by Nj C n the set containing 
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indices of all the external edges connected to the j-th one by an inner edge, i.e. 

Nj := {kEm: S jk / 0}U{fc6m: (3/ > m + l)(T jt + A T kl ^ 0)} 

U {k > m + 1 : Tjk ^ 0} for j G rh 

Nj := {kern: T kj ^ 0} for j > m + 1 

The definition of the set Nj has two simple consequences, namely 

k e Nj j e N k and j>m + l=>Nj<Zth. 

We employ the following symbols for wave function components on the edges: those 
on the j-th external one is denoted by fj, while the wave function on the connecting 
segments is denoted f{j,k) on the interval between V{j±y and Vj and f( k ,j) on the other 
half of the segment; the conventions about parametrisation of the intervals have been 
specified above. 

Next we shall write explicitly the coupling conditions involved in the above described 
scheme, first without the vector potentials; for simplicity we will often refrain from 
indicating the dependence of the parameters Wj(d), w^j ik y(d) on the distance d. The 5 
interaction at the segment connecting the j-th and k-th outer edge (present for j,k G n 
such that k G Nj) is expressed through the conditions 

/a,*) (o) = = : W°) > + = w {j,k}fim(°) > 

while the 5 coupling at the endpoint of the j-th external edge, j G h, means 
fM = fm(d) f o r all k G Nj, fj(0) - £ f m (d-) = Wjfj(0) . 

k&Nj 

It is not difficult to modify these conditions to include the vector potentials using a 
simple gauge transformation |CET l0j: the continuity requirement is preserved, while 
the coupling parameter changes from Wj(d) to Wj(d) + iJ2 k eNj A(j )k )(d); in other words, 
the impact of the added potentials results into the phase shifts dA(j tk ){d) and dAi kt j\(d), 
respectively, on the appropriate parts of the connecting segments. 

Using the above conditions one can find suitable candidates for Wj(d), W{j )k }(d), and 
Au 7 k)(d) by inserting the boundary values written as 

f m (d) = e idA ^ {f m (0) + df [hk) (Q)) + 0(d 2 ) and 
f' m (d) = e idA ^ f' m (0) + 0{d) 

for any j,k G h and fixing the ^-dependence in such a way that the limit d — > 
yields (12.21) . The procedure is demanding and described in detail in CETlOj . we will 
mention just its results. As for Arj )k Ud), we have the relations 

A Qk) (d)=U^ TgTjk , *X»T sk >0, 
U ' k)K ' \ i (arg T jh - tt) if Re T jk < 1 ; 

for all j G rh, k G Nj \ rh, while for j G rh and k G Nj fl rh we put 

= W r A/c _ , A t ^ _ j ( 2 - 3b ) 



2</ 



depending similarly on whether Re (dSj k + ^™ =m+1 Tj{T k i) is nonnegative or not. Con- 
cerning W{j )k y(d), we require that 

w {j,k}( d ) = \{^ 2+ (7^)") V j G m, k e Nj\rh . (2.3c) 
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and 

1 

2 + d ■ W{j, k } 



-ld-S jk + T ji T ki ) v 3 em, keNjnm, (2.3d) 



l=m+l 



where the bracket symbol on the right hand side is defined as (c) := |c| for Rec > and 
(c) := —\c\ for Rec < 0. Finally, the expressions for w k are given by 

w k (d)= LziM±pkgk) Vfc>m + 1, (2.3e) 

Lb 

and 

I AT I m 1 n 1 " 

«>i(<o = % + d E c 1 + w ( 2 - 3f ) 

fc=l Z=m+1 i=m+l 

if j G m and k E Nj f] rh. 

Remark 2.2. For our later considerations it is crucial to know precisely the dependence 
of the magnetic and electric potentials A^) = ^0',fe)(^) an d w {j,fc} = w {j,k}(d) on the 
internal length d. We have = ff(d~ 1 ) and w^^d) = &{d~ 1 ) if fc G iV,- \ m. If 

G iVj fl m, then we have to distinguish two cases. If 

■n 

T ^ ± °> ( 2 - 4 ) 

Z=m+1 

then we again have w^^id) = G(dr x \ Otherwise, we collect another power of d" 1 and 
obtain wijk}{d) = @{dr 2 \ We are not aware of any meaning of (j2.4p in terms of the 
original vertex coupling or equivalent characterisations. 

The choice of the parameters has been guided by formal considerations but it opens 
way to prove the convergence of the corresponding operators. Let us denote the Laplacian 
on the star graph T(0) with the coupling (I2.2p in the vertex as H stav , while // approx w ill 
stand for the operators of the described approximating family; the symbols i? approx (z) 
and i?^ pprox (z) will denote respectively the resolvents of those operators at the energy 
z outside the spectrum. We have to keep in mind that they act on different spaces: 
R st&r {z) maps L 2 (r(0)) onto dom# star , while the domain of Rf piox {z) is L 2 (r 5 ' T (d)), 
where T s ' T (d) = T(0) U Tf^(d) and where T i ^(d) is the graph of connecting (inner) 
edges of length 2d described above. In order to compare the resolvents, we identify thus 
R sta,: (z) with the orthogonal sum 

RT(z) ■■= R st&v (z) ©0 (2.5) 

adding the zero operator acting on L 2 (T^(d)). Then both operators act on the same 
space and one can estimate their difference; using explicit forms of the corresponding 
resolvent kernels one can check in a straightforward but rather tedious way the relation 

\\R^(z)-Rf pmx (z)y 2 = 0(Vd) as d -> 0+ 

for the Hilbert-Schmidt norm, see jCET IO]. With the identification ([23]) in mind we can 
then state the indicated approximation result. 

Theorem 2.3. Let Wj(d), j G n, wij >k }{d), j G h, k G Nj and A^ ,k \d) depend on the 
length d according to fl2.3ajl - f!2.3fj) . Then the family 7f approx converges to if star in the 
norm-resolvent sense as d — > 0+. 



We present some examples of vertex coupling approximations in Section 15. 2[ 
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3. Approximation by Schrodinger operators on manifolds 

Now we pass to the second step and show how the intermediate quantum graph con- 
structed in Section [2] with 5 couplings and vector potentials can be approximated by 
scaled magnetic Schrodinger operators on manifolds. For the sake of simplicity, we con- 
sider first an approximation using abstract manifolds without boundary, and discuss the 
case of a graph embedded in R u subsequently in Section 15.11 To set up the approxima- 
tion scheme, it is convenient to work with appropriate quadratic forms instead of the 
associated operators. 

3.1. The spaces and quadratic form on the graph level. We start with the defi- 
nition of the Hilbert space and quadratic form on the intermediate graph T = T s,T (d), 
where d £ (0, 1] denotes the approximation parameter of the previous section. It is con- 
venient to modify slightly the convention (jrv]) concerning the internal edges e = {j, k}; 
from now on we shall consider each of them as a single edge with the 5 interaction in the 
middle (i.e. at vy^}) an d identify this edge with the interval [— d, d], oriented in such a 
way that the parameter increases from j to k if j < k. Concerning the vector potential, 
we set A e := Au^ = —Aujy For the sake of brevity, we use the symbols A = (A e ) e , 
w = (w e ,w v ) etV for the collections of magnetic potentials and S interaction strengths, 
respectively. We will also often suppress in the sequel the dependence of the quantities 
on d, A, and w. With each outer edge e £ n = {1, . . . , n}, we associate I e := [0, oo), and 

for each inner edge e e g) = { {] ,k } 1 1 < j < k < „ }, we set /. = I M = As 

the Hilbert and Sobolev spaces on a fixed edge needed in our approximation we set 

:= L 2 (I e ) and J^ 1 := H 1 ^), 

where L 2 (J) and H 1 (7) denote as usual the space of square integrable functions and of 
once weakly differentiable and square integrable functions on the interval J, respectively. 

For all the quadratic forms defined below, the domains consist of elements of M^. With 
the described parametrisation of an inner edge e = {j, k} with i < k the corresponding 
quadratic form is 

Ufe) ■= f d \f'e(s) + iAJ e ( S )\ 2 ds + W e \f e (0)\ 2 - 
J-d 

This form corresponds to the Laplacian on the edge with the magnetic potential A e and 
the S interaction at the point s = 0. It is convenient to introduce also a quadratic form 
which includes the effect of the 5 interactions at the edge endpoints, namely 

Ufe) := Ufe) + ' |/eM)f + 777~T ' \fe{d)\\ 

\^j\ 

On an outer edge, we simply set 

Ufe)--= / \f'e(s)\ 2 d S . 
JO 

The full Hilbert and Sobolev spaces are 

:= J% and J^ 1 := n C(r), 

e e 

where the sum runs over all the inner and outer edges. More explicitly, the Sobolev 
space J$? 1 consists of all functions in H 1 (/ e ) on each edge, which are continuous on T, 
i.e. which have a common value 

{/ e (0), if e = j is an outer edge, 
fe(-d), if e = {j, k} ~ v = Vj is an inner edge, j < k, 
f e (d), if e = {j, k} ~ v = v k is an inner edge, 
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for all edges e ~ v, i.e. adjacent with v. 

The quadratic form on the intermediate graph T(d) is given by 

f)(/) :=5>(A) 

e 

for / = (f e )e £ JF 1 ] the corresponding operator is the one described in Section [2] with 
S interactions of strength Wj at vertex Vj and of strength w e in the middle of the inner 
edge e = {j, k}, as well as vector potential Aq^) supported by this edge. 

For comparison reasons, we also need the free quadratic form, without both the mag- 
netic potentials and the 6 interactions, which is given by 

Mfe):= I \fe{s)\ 2 ds and := 

with the same domains as f) e and f), respectively. It is easy to see that is a closed 
quadratic form, i.e. that domfl = with the norm given by ||/||^i := ||/|| 2 + 'O(f) is 
complete, and therefore itself a Hilbert space. The operator corresponding to t> is the 
free Laplacian on T(d), often also called Kirchhoff Laplacian on the graph. 

Proposition 3.1. 

(i) The quadratic form f) is relatively form-bounded with respect to d with relative 
bound zero. More precisely, for any r\ > there is a constant C v > depending 
only on r\, d, A := max e \A e \, andw := 3 max e ^{|w e |, \w v \} such that 

\i)(f)-Hf)\<vHf) + cjf\\ 2 . 

In particular, f) is also a closed form. 

(ii) Wehaved(f)<2(t)(f) + C 1/2 \\f\\ 2 ). 

Proof. §j§ On the interval [—d, d] we have the following standard estimate 

|/Wr<a||/T + ^II/H 2 (3-1) 

1 1 a 

for all s G [—d, d], < a < d, and / G H 1 (— d, d). Moreover, for any rj > and a, b G K 
we have 

^•^-^^(a + ^^a + ^.^+^ + I).^. (3.2) 
In particular, for an inner edge e = {j, k} we have 

Uf) - Uf) = \\f + iA/|| 2 - ll/'f + w e \f(0)\ 2 + ^\f(-d)\ 2 + ^\f(d)\ 2 

on [— d, d] using (13.1 p with s G {— d, 0, d} and the upper estimate in (13.21) with rj/2 instead 
of rj, where 

— ii, \ W j\ \ W k\ 

W e := \w e \ + r-r-r + 



\Nj\ \N k 



a := min<j — — , d }> (3-3) 



Choosing 

a ■= min <! 

we can estimate the coefficient of e (/) — ||/'|| 2 by r\. 
For the opposite inequality, we have 

■2\A e \ 2 2w { 



Uf)- US) ^(i-T^ + ^) ii/'f + (^+ „ 
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using now the lower estimate in ( 13. 2 p with rj/2. In particular, with a as in (I3.3P and with 
1 — (1 + r//2) _1 < 77/2 we can again estimate the coefficient of e (/) = ||/'|| 2 by rj. As 
constant C nfi on each edge, we can therefore choose 

^ e :=(l + ^)|^| 2 + max{^,^}. (3.4) 

Summing up all contributions for each edge, we can choose C v := max e C^ e , and this 
constant depends only on rj, d, A and w. 
(pi]) follows with rj = 1/2. In particular, 

C 1/2 = C 1/2 (gU,u;) = ^(Z 2 ) + 0{w 2 ) + ^(|)- (3.5) 

□ 

3.2. The spaces and quadratic form on the manifold level. We now define the 
manifold model as in [EP09] . For a given e G (0, d] we associate a connected (m + 1)- 
dimensional manifold X £ to the graph r(d) as follows: To the edge e and the vertex v 
we associate the Riemannian manifolds 

X £fi := I e x eY; and X £>v := eX„, (3.6) 

respectively, where eY e is a manifold Y e of dimension m > (called transverse manifold) 
equipped with the metric h £>e := e 2 h e . More precisely, the so-called edge neighbour- 
hood X £ie and the vertex neighbourhood sX £tV carry the metrics g £fi = d 2 s + e 2 h e and 
9e,v = ^ 2 5'd, where h e and ^ are e-independent metrics on Y e and X v , respectively. We 
assume that for each edge e adjacent to v, the vertex neighbourhood X £>v has a boundary 
component d e X £>v = ed e X v isometric to the scaled transverse manifold eY e . Fixing such 
an isometry and assuming that X £jV has product structure near each of the boundary 
components d e X £jV , we identify the boundary component d v X £fi = {0} x eY e of the edge 
neighbourhood X £i£ with d e X £>v . 

For simplicity, we assume here that the transversal manifold Y e has no boundary and 
that its volume is normalised, i.e. vol m Y~ e = 1. 

On a Riemannian manifold X, we denote by L 2 (Jf) the Hilbert space of square inte- 
grable functions on X with respect to the natural measure induced by the Riemannian 
metric. Moreover, we denote by H 1 (X) the completion of the space of smooth functions 
with compact support (not necessarily vanishing on the boundary of X) with respect 
to the norm given by ||w||hi(x) := IMIl 2 (x) + II^ m IIl 2 (x)' where du denotes the exterior 
derivative of u on X. 

We set 

J^ >e := L 2 (ie, Jf £l e), X~ £ , e :=\- 2 (eY e ) and M^ tV := l 2 (X £>v ). 

We will often identify an L 2 -function u on X £j6 with the vector-valued function I e — > J^ >e , 
s 1 — y u(s) := u(s, •). 

For each inner edge, we set 

fk> e ) := J (||w' e ( s ) + iA e Ue(s)\\ 2 + t £>e (u e (s))) ds + ^ J ||u e (s)|| 2 ds, 
where u' e denotes the derivative with respect to the longitudinal variable s and where 

teM := \\ d YM\l 2 {eY e y 

Here, dy e y9 is the exterior derivative on the manifold Y e . For each outer edge we set 

POO 

M«e) ; = / (\K(s)\\*r ee + t e>e (u e (s)))ds 
Jo 
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In both cases, u e G J^f e \ = H 1 (X £je ). On a vertex neighbourhood, we set 

2 i Wv II II 2 



evo\x v ]]Uv]]L ^y 

The total Hilbert spaces here are 

^:=0^, e ©0^ and J^ 1 := H\X £ ), (3.7) 

e v 

where the sum runs over all inner and outer edges. Now, the quadratic form on the 
manifold X £ is given by 



e v 



for u G Jrf? 1 with the obvious notation u e := u\ Xge and u v := u\ Xev . The corresponding 
operator is a magnetic Schrodinger operator on X e with (constant) potential w e / (2s) on 
[— £, e] x eY e in the middle of an edge neighbourhood and w v /(evo\X v ) on each vertex 
neighbourhood. For the use of non-constant potentials we refer to [EP09J. 

For comparison reasons, we also need the free quadratic form (i.e. without magnetic 
and electric potentials), given by D £te {u e ) = ||dw e || L 2 (jvr e e )> ^£,f(u v ) = \\<iu v \\^ Xe ^ and 

e v 

with the same domains as for f) £ie , \) £jV and f) e . Since we define J)^ 1 = H^X^) as the 
completion of smooth functions with compact support with respect to the norm ||w||^i := 
9 £ (m) + ||m|| 2 , the quadratic form £ is closed. The operator corresponding to X) is the 
Laplacian on X £ . 

Proposition 3.2. 

(i) The quadratic form f) e is relatively form-bounded with respect to d £ with relative 
bound zero. More precisely, for any r\ > there is a constant C v > C v > 
depending only on r\, d, A : = max e \A e \, w := 3 max e t; {|w e |, \w v \} and X v such 
that 

\t) e (u) - D £ (u)\ < r]d £ (u) + CJu\\ 2 (3.8) 

for all < e < Eq, where e := t)c(v)/\w v \ and where c(v) is a constant depending 
only on X v . In particular, f) £ is also a closed form, 
(li) We have D £ (u) < 2(fj e (u) + Ci/ 2 \\u\\ 2 ) . 

Proof. The proof is very similar to the one of Proposition 13.11 For fli]), we have the 
following vector- valued version of (13. ip . namely, 

l|2 2 
IK^)^ < a|Kll^j, e + -||«e||jft ia (3-9) 

for all s G [—d, d], < a < d and u G H 1 (X e e ). In particular, for an inner edge e = {j, k} 
we have 

|f) e ,e(We) - *MO| < V\K\\%, e + C ri,e\\ U e\\% £<e 

with C Vte as in (13.41) . 

On a vertex neighbourhood, we have 



I t ( \ ( \ I I W I II 1 1 2 



W„| 
£ Vol X„ 
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for < a < d using |EP09[ Lem. 2.9], where c wo \{y) := vol X v j vol m dX v and C{y) is 
another constant depending only on X v , see jEP09] for details. Setting 

a := mm{d,r]vo\X v /(Ac vo i(v)\w v \} and e :=mm- — - " , 

v \w v \C{v) 

and summing up all contributions, we can choose C v > such that (13. 8 p holds for all 
< e < e with 

C v = C v (d,A,w) = 0(A 2 (l + -)) + ) + (3.10) 

and the error term depend additionally only on X v . The remaining assertion fln|) follows 
as before. □ 

4. CONVERGENCE OF THE OPERATORS 

4.1. Norm convergence of operators and forms acting in different Hilbert 
spaces. Let us briefly review the concept of norm convergence of operators acting in 
different Hilbert spaces introduced first in |P06t App.]. A general spectral theory for 
quasi-unitary equivalent operators is developed in a more elaborated version in [P12, 
Ch. 4], see also [EP09] . 

Let and 1 be Hilbert spaces such that Jif is a dense subspace of Jff with 
||/||.# < ll/lljr 1 an d similarly for C 3^. Let \) and f) be closed, quadratic forms, 
semi-bounded from below with domain and J^ 1 , respectively. 

Let 5 > 0. We say that f) and f) are 5-quasi-unitarily equivalent if there are so-called 
identification operators 

J: 3¥ — > j¥, J 1 : Jf 1 — yJT 1 and J' 1 : JF 1 — JT 1 , 

such that these operators are 5-quasi unitary, i.e. 

II Jf ~ J'ff < S 2 \\f\\%„ \\J*u - J' l u\\ 2 < (4.1a) 
|| rjf _ ff < ^H/H^, || JJ*n - M || 2 < ^IMI^, (4.1b) 

<5|H|^||/||^ (4.ic) 

for / and u in the appropriate spaces. The attribute 5 -quasi-unitary refers to the fact 
that we have a quantitative generalisation of unitary operators. In particular, if 5 = 0, 
then a 5-quasi-unitary operator is just unitary. 

On the operator level, we have the following definition: Denote by H and H the (self- 
adjoint) operators associated to () and f). We say that H and H are 5-quasi-unitarily 
equivalent (see again Footnote [T|) if there is an identification operator J : Ji? — > Jf? 
such that 

IKid-JV)^!! < 5, || (id -JJ*)^ || < 5 and || JR+ - j\\ < 5, (4.2) 

where ||-|| denotes the operator norm, and where := (H =p i)" 1 and R^ 1 := (H =p i) _1 
denote the resolvents, respectively. The resolvent estimates are supposed to hold for 
both signs. 

We have the following relation between the quasi unitary equivalence for forms and 
operators. For convenience of the reader, we give a short proof of the first assertion here. 
The remaining assertions follow from the abstract theory developed in [P06t App. A] 
and [ELI Ch. 4]. 



1 We warn the reader that in [PI 2) the notion "<5-quasi-unitary equivalent" is defined in a slightly 
more general way (allowing e.g. a second identification operator J' : — > Jf? such that || J* — J'\\ < 5 
to cover some more general situations). This should not cause any confusion here. 
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Theorem 4.1. Let 5 > and C > 1. Assume that I) and f) are 5-quasi-unitarily 
equivalent closed quadratic forms such that 

11/11^ <2(f)(/) + C||/f) and \\u\\ 2 ~ 1 <2(l ) (u)+C\\u\\ 2 ) 

for all f G and u G M 31 . Then the following assertions hold: 

(i) The associated operators H and H are (l2Cd~)-quasi-unitarily equivalent. 

(ii) There is a universal constant c(z) > depending only on z such that 

\\J(H - z)- 1 -(H- z)- l J\\ < c(z)C5, (4.3a) 
\\J{H - z)- 1 J* - (H £ - z)' 1 \\ < c{z)C5 (4.3b) 

for z G C \ R. Moreover, we can replace the function <p(\) = (A — z)" 1 in 
(p(H) = (H — z)~ l etc. by any measurable, bounded function converging to a 
constant as A — >• oo and being continuous in a neighbourhood of a(H). 

(Hi) Assume that H = H £ is 5 e -unitarily equivalent with H , where 5 e — > 0, then the 
spectrum of H £ converges to the spectrum of H in the sense that if \ £ G o~(H £ ) 
and X £ —> X, then A G o~(H), and if A G cr(H), then there exists (\ £ ) £ such that 
X £ G o-(H £ ) and X £ — 0. The same is true for the essential spectrum. 

(iv) Assume as before that H = H £ is 5 £ -unitarily equivalent with H , where 5 £ — > 0, 
then for any X G crdisc(H) there exists a family {X £ } £ with X £ G (Xdisc(-ffe) such 
that X £ — > X as e — > 0. Moreover, the multiplicity is preserved. If X is a simple 
eigenvalue with normalised eigenfunction if, then for e small enough there exists 
a family of simple normalised eigenf unctions {<p £ } £ of H £ such that 

\\J<£ ~ <£> e ||L 2 (Xe) -> 

holds as e — )■ 0. 
Proof. ^ From our assumption, we have 

ll/ll,^ < mi) + C\\f\\ 2 ) = 2|f,(/) + ||/|| 2 | + 2(C - 1)||/|| 2 . 
Moreover, the first term can be estimated as 

|t)(/) + ii/irr<2(f)(/) 2 +ii/ii 4 ) 

= 2|[)(/)-i||/|| 2 ||h(/) + i||/|| 2 | 

= 2|((ir=Fi)/,/>||</ l (ir=Fi)/>| 
<2||/|| 2 ||(iJ T i)/|| 2 <2||(iJ T i)/|| 4 

using || (if =p i) _1 || < 1 at the last step. In particular, we have 

II/U3,* < (2y/2 + 2C- 2)\\(H T i)/|| 2 < AC\\(H T i)/|| 2 (4.4) 

since 2\/2 - 2 < 2 < 2C. Similarly, we can show the same estimate for u, and we have 

\\f\Ui <2VC\\(H T i)f\\ and < 2VC\\(H T i)w|| ■ (4.5) 

Therefore, we conclude 

II/- JV/H < 8\\f\\jti < 2VC8\\(H =Fi)/|| 

by (14.1bj) . and in particular, ||(id — J*J)R ± \\ < 2\J~Cb~. The second norm estimate in ( 14. 2 p 
follows similarly. 
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For the last norm estimate of the quasi- unitary equivalence of the operators in ( I4.2|) . 
set / := R g G domH and u := R T v G domif. Then we have 

((JR ± - R±J)g, v) = (Jf, v) - (g, J*u) 

= ((J - J 1 )/, v) + ((J 1 /, (H ± i)u) — ((H =F i)f, J' l u)) 

+ (g, (J n - r)u) 
= ((J - J 1 )/, v) + (f)( J 1 /, u) - ()(/, J n u)) + ( 5 , (J' 1 - J» 

T i(((ji-j)f, u } + (f,(r-j' 1 ) u }), 

and therefore 

(((J^-^J)^,^)! < (2VC + AC + 3 ■ 2VC)S\\g\\\\v\\ < 12C6\\g\\ \\v\\ (4.6) 



using ( l4~Tj) and (145]) . 

Once we have the estimates of the quasi- unitary equivalence in (14. 2p . the remaining 
assertions follow as in |P06t App. A] or |P12[ Ch. 4]. □ 

We remark that the convergence of higher-dimensional eigenspaces is also valid, how- 
ever, it requires some technicalities which we skip here. 

Remark 4.2. Note that we only obtain the quasi-unitary equivalence of the operators 
with a factor C and not y/C. This is due to the fact that from (14. left , we collect 
two factors 2\fC for the estimates H-R^gH^i < 2\/C'||^||jf and \\R T v < 2y/C\\v\\jp 
in (1461) . 

4.2. Quasi-unitary equivalence between the graph and manifold forms. We 

now apply the abstract results of the previous section to our problem where 

J? := l 2 {T s ' T {d)), je 1 := H^r^d)), JF:= L 2 (X £ ), JF 1 := H\X £ ). (4.7) 

We start with the definition of the identification operator on an edge. Let 

J e : Jf e = L 2 (J e ) — > J^ £ ,e = L 2 (A £>e ) be given by J e f e := f e ® l £je , 

where l £>e is the (constant) eigenfunction of Y e associated to the lowest (zero) eigenvalue 
equal to e~ m ^ 2 . Since we assumed vol Y" e = 1, the eigenfunction is normalised. Its adjoint 
acts as transverse averaging, 

(J* e u e )(s) = (u e (s), l £l eU, e = e m/2 / u e {s,y e ) dy e . 

Before defining the global identification operator, we need the following result: 
Lemma 4.3. For 0<d<l ; 0<£<l and f,g G H 1 ([— d, d\) we have 

1 ' £ f(s)g(s)ds-f(0)g(0) < 2(e/d)^\\f\\ H1 \\g\\ Hi . (4.8) 



i/(s)r<~ii/ir H1 (4.9) 



2e 

Proof. Note first that 



for s G [—d, d] by (13. ip since d G (0, 1] by assumption. From f(s) — /(0) = J* f'(t) dt we 
conclude 

i/( S )-/(o)r<| S |iini 2 . (4.io) 
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Now, the left-hand side of (14. 8p can be estimated by 
Y e f\m-m\\9(s)\ds+ 1 -^- Jjg(s)-g(0)\ds 

^^ll/llyfNlHxv^+iy^ll/llH^lb'll 

using fl4.9p - fl4.10p together with Cauchy-Schwarz inequality, from where the desired es- 
timate follows. □ 

We can now compare the two contributions of the quadratic forms on an internal edge, 
including the potential in the middle of this edge. We could consider this inner point as 
a vertex, too, and use the arguments for vertex neighbourhoods as in [EP09J. Since this 
vertex has degree two only, we give a direct (and simpler) proof here: 

Lemma 4.4. We have 

\t)e,e(Jefe,U e ) - t) e (f e , J*U e )\ < 2 1 W e \ {s / d) 1/2 || / 1| H i(r) \\u || H i(X e ) 

for all f G Jrf? 1 = H 1 ^), ueif i = rfpf £ ) ; < e < 1 and < d < 1. 
Proof. We have 

t)s,e( J efe, U e ) ~ t)e(fe, J* e U e ) 

\(f' e ® l £ , e + \A e f e ® l eje ){ S ),U e {s))jt e , e ~ (f'e(s) + i4/ e («))(a e (s),l E ,e>^) d S 

+ W e (J (fe(s)l £>e , U e { S ))^ e ds - / e (0) (u e (0), l e , e )jr e , e ) ■ 

Note that in the first integral the term with the derivatives and the magnetic potential 
contributions respectively cancel. Moreover, the expression contains no contribution 
from the transversal (sesquilinear) form fi £>e since t £ie (l £je , </?) = for any y> G L 2 (eY" e ). 
The remaining (electric) potential term can be estimated by Lemma I4T31 with / = f e and 

g{s) = (u e (s), l e , e )jr ei J. □ 

As the global identification operator we define J : — > Jff by 



with respect to the decomposition (13. 7p . In order to relate the Sobolev spaces of order 
one we correct the error made at the vertex neighbourhood by fixing the function to be 
constant there. Namely, we define J 1 : Jrf? 1 — > Jif 1 by 

J l f:=@Jefe@e- m/2 0/(v)l„, 

e v 

where t v is the constant function on X v with value 1. Since / is continuous on the graph, 
J f is continuous along the vertex and edge neighbourhood boundary, and therefore maps 
into the Sobolev space J^ 1 = H 1 (X e ). 

For the operator J' 1 : ^C 1 — > Jff 1 , we have to modify J* in such a way that the first 
order spaces are respected, namely we set 

(J?u)(s) := (JX)( a ) + X-{s)e m ' 2 (j Vj n - (J> e )M)) 

+ X + (s)e m/2 (f u-{Jtu e ){d)) 



APPROXIMATION OF VERTEX COUPLINGS BY SCHRODINGER OPERATORS 



15 



on an inner edge e = {j, k}, j < k, where x± are smooth functions with x±(±cf) = 1, 
|x±l — 2/d and X±( s ) = f° r ^ s — 0- Moreover, 

f* u := ^, { ^ U) := ^x-J Xv UvdXv 

is the average of a function w on the (unsealed) vertex neighbourhood X v . 
On an outer edge e = j we set 

:= G>X)00 + x(^ m/2 (.f> - (J>e)(0)) 

where \ is a smooth function with x(0) = 1, |y/| < 2 and x( s ) = f° r s > 1. Note that 
J' 1 differs from J*f only by a correction near the vertices. Since (J' 1 w) e (f) = E m l 2 j- v u 
independently of e ~ v, the function J n u is indeed continuous, and therefore an element 
ofH 1 ^). 

Now we can make a claim which represents the main technical ingredient in the analysis 
of the two quadratic forms: 

Proposition 4.5. Let < d < 1, then the quadratic forms f) e and f) are 5 £ - quasi-unitary 
equivalent, where 5 £ depends on e, d and w := 3 max e ^{|iy e |, \w v \} as follows 

«, = ,((«)>+!))+,(!£). 
Moreover, the error depends additionally only on X v and Y e . 

Proof. The argument is similar as in the (simpler) situation of [EP09 , Prop. 3.2] using the 
identification operators just defined. A new feature here is that we employ additionally 
Lemma 14.41 to compare the form contribution on the internal edge neighbourhood and 
its counterpart on the metric graph. Note also that in the present situation we have 
additionally the magnetic potential and slightly different constants than in [EP09] . □ 

In particular, by a clever choice of the e-dependency of the parameter d, we are able 
to make the following conclusion: 

Corollary 4.6. Assume that w e , w v and A e are chosen as in (I2.3p . then w = €?(d~ 2 ) 
and A = G{d~ x \ If in addition, d = e a with < a < 1/5, then f) e and f) are S e -quasi 
unitarily equivalent for all < e < e\, where 5 £ = ff{e^ l ^ a ^ 2 ), and where E\ > is a 
constant. 

Finally, if < a < 1/13, then the associated operators H e and H are S £ -quasi unitarily 
equivalent with 5 E = ^(£ (1 " 13q)/2 ). 



Proof. The quasi- unitary equivalence of the quadratic forms follows from Proposition I4.5[ 
as well as the estimate on S e . Moreover, Eq = Eq(e) as given in Proposition 13 . 21 is generally 
of order 0{\jw) = 0(E 2a ), i.e. e < c£ 2a . In particular, we can choose, E\ = c 1 ^ 1 ^ 2 ^. 
For the last assertion, note that the constants Cyi and Ci/ 2 of Propositions 13.11 

and O fulfil Ci /2 = @(E~ Aa ) and C 1/2 = ^(e" 4a ), cf. §£5$ and f l3~TU|) . since the term 
w 2 = ff{E^ i ) is dominant. The result now follows from Theorem 14.11 fli]) with C : = 

max{C 1/2 , C 1/2 }, and therefore we have \ = 12C5 £ = ^( e -4»f(i-5«)/2) = ^( e (i-i3a)/2^ 

□ 

Finally we are in position to put the two convergence steps together; and to state and 
prove the main result of this article: 

Theorem 4.7. Assume that T(0) is a star graph with vertex condition parametrised by 
matrices S and T as in Section^ and let < a < 1/13. Then there is a Schrodinger 
operator H e on an approximating manifold X £ as constructed in Section VS.'lA such that 

\\JRf aX (z)J* -R e (z)\\ = ^>( £ niin{l-13a,a}/2) 

for z e C \ R, where R £ {z) = (H £ -z)~ l . 
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Proof. The result is an immediate consequence of Corollary I4.6[ Theorem 14.11 (jn|) and 
Theorem 12.31 □ 

Remark 4.8. The error term in the theorem depends only on z and the building block 
manifolds X v at the vertices and the transversal manifolds Y e on the edges. If a = 1/14, 
we obtain the error estimate ^(e 1 ^ 28 ) which is the maximal value the function a \— > 
min{l — 13a, a}/2 can achieve. 

The error estimate we obtain here is of the same type that we obtained in [EP09| Sec. 4] 
when we approximated the S's interaction despite the fact that the present approximation 
of this particular coupling is different, cf. Section l5\2l below. 

If the condition (12.41) mentioned in Remark 12.21 is fulfilled for all j, k we obtain a 
slightly better estimate. In this case, we have w = G'{d~ 1 ) instead of &[d~ 2 \ and f) e 
and f) are <5 £ -quasi, where 5 £ = ff{e^ l ~^ a ^ 2 ). Moreover, the associated operators H e and 
H are 5 e -quasi unitarily equivalent with 5 £ = ^(e^ 1-7 "^ 2 ). However, both assumptions 
made about a, namely < a < 1/13 and < a < 1/7, are for sure not optimal. 

There is an obvious extension to the above convergence result for quantum graphs T 
with more than one vertex. For quantum graphs with finitely many vertices, the con- 
vergence result holds without changes, and for infinitely many vertices, some uniformity 
conditions are needed. Such questions are discussed in detail in [P06j and |P12j . 

Remark 4.9. One may ask whether one can reformulate the "quasi-unitary equivalence" 
for the present situation using 

J, L 2 (r(o)) — ► i 2 (r s ' T (d)) = L 2 (r(o)) © L 2 (rf nt T (rf)), If = f © o, 

in which case R s ^ av (z) = JR star (z)J* by (12. 5p and the resolvent convergence of Theo- 
rem [2]3] can be stated as 

\\JR stax (z)J* - i?r rOX (^)IU(L 2 (r^ W) ) = *(d 1/a ) (4-11) 

for d — > 0. In fact, we are interested primarily in spectral consequences of such a 
reformulation which can be demonstrated in a more direct way. To this end, note that 
eq. (14. lip is just (14.3bj) of Theorem 14.11 without the constant C . Moreover, from |P12l 



Thm. 4.2.9-10] one can conclude that (I4.3al) is valid for more general ip than (p(X) 
(A — z)~ x , see Theorem 14.11 ([II]). Using arguments analogous to those in |P12l Sec. 4.2- 
4.3], we can deduce from (14.111) that (14.3b|) also holds for such ip. Consequently, the 
spectral convergence stated in Theorem 14.11 (lull) and flivl) also holds in this situation. 

5. Examples 

5.1. Embedded graphs and graph neighbourhoods. Consider the situation when 
the graph is embedded in M. u , v > 2. This may be a restriction to the vertex coupling 
if v — 2 and the vertex degree exceeds three; recall that the edges of the internal graph 
defined in flmj) of Section [T] are supposed to be non- intersecting. For v > 3 this difficulty 
can be avoided in the edges are properly curved. At the same time, irrespective of d the 
lengths of the edge parts of the manifold change as e — > by an amount given by the 
size of the vertex neighbourhoods. Let us point out briefly that for such embedded "fat 
graphs" curved and shortened edges lead to a small error in the approximation only. 

Consider first the length change. In our case, the difference of the original edge length 
and the one of an embedded edge is of order d — e = d{l — e/d) = d(l — e l ~ a ). We 
have shown in |EP09| Lem. 2.7] that this leads to an additional error of order ff{e 1 ^ a )\ 
expressed again in terms of quasi-unitary operators. 

Furthermore, if we allow curved edges in the case of a graph embedded in W, we still 
arrive at the same limit operator. The error is of order ^(£ 1-Q ) (see |P12l Sec. 6.7 and 
Prop. 4.5.6] for details; the factor e comes from the shrinking rate, the factor e~ a from 
the curvature term of the embedded curve in dimension v = 2; the length shrinks by 
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d = s a so its curvature is of order e~ a ). Similar arguments apply for v > 3. In particular, 
combining the effect of shortening of edges and curved edges, and using the transitivity 
of quasi-unitary equivalence ( |P121 Prop. 4.2.8]) we arrive at an error estimate which is 
not worse than the one in Theorem 14.71 

5.2. Special vertex couplings and approximation by Schrodinger operators. 

While the approximation described in Theorem 14.71 cover any self-adjoint coupling, for 
some of them we have better alternatives. This concerns, in particular, the 5 coupling 
where a simple scaled potential does a better job as explained in [EP09] . On the other 
hand, for couplings with functions discontinuous at the vertex we do not have many 
alternatives. 

It is illustrative to compare the approximation of the 5' s coupling obtained from 
the graph-level approximations described in Section [2] with the one from [CE04| used 
in |EP09] for the approximation by Schrodinger operators. Recall that a 5' s coupling of 
strength (3 in a vertex of degree n edges characterised by the condition 

^/(o) - /'(0) = 0, 

where J is the n x n matrix with all entries one. In other words, the respective ST- 
parametrisation from Proposition 12.11 is given by m = n, S = ft' 1 J and T = 0, and the 
strengths of the 6 potentials required to approximate S' s according to Theorem 12.31 are 

P 2 , 2-n n-l 
W{ ^ = -J>-d and W > = — d~- 

In particular, all inner edges are present. If n — 3, for instance, we employ a small 
triangle graph of length scale d = e a attaching the "external" edges to its vertices (as 
sketched in Figured]). The corresponding Schrodinger operator has a potential of order 
_ £ -«-i near v . anc [ f order — fie~ 2a ~ x at the midpoint of each edge {jk}; for simplicity 
the potentials can be chosen piecewise constant. 

The approximation used in [EP09J is different. Here we keep the original star graph, 
but introduce additional 5-couplings on each edge at distance d = e a of the central vertex. 
The strength of the coupling at the central vertex is —(3/d 2 , hence the Schrodinger 
potential there is of order — (5e~ 2a ~ l . The strength of the coupling at the additional 
vertices is —l/d, hence the Schrodinger potential is of order —f3e~ a ~ l . One sees that the 
approximation graph topology is different but the S strengths in the two cases differ only 
in lower order termg^ with respect to the length scale d = e a . 

Let us finally remark that S' s is not the only example of interest; our method makes 
it possible to approximate other couplings of potential importance such as the scale- 
invariant ones analysed recently in |( 'KT1 lj. 
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